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Abstract. We aim at an explicit characterization of the renormalized Hamiltonian after dec- 
imation transformation of a one-dimensional Ising-type Hamiltonian with a nearest-neighbor 
interaction and a magnetic field term. To facilitate a deeper understanding of the decimation 
effect, we translate the renormalization flow on the Ising Hamiltonian into a flow on the asso- 
ciated Markov chains through the Markov-Gibbs equivalence. Two different methods are used 
to verify the well-known conjecture that the eigenvalues of the linearization of this renormal- 
ization transformation about the fixed point bear important information about all six of the 
critical exponents. This illustrates the universality property of the renormalization group map 
in this case. 
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1. Introduction 

The discovery of the equivalence of Markov random fields and Gibbs random fields was a 
major breakthrough in the interchange of ideas between probability and physics. A Markov 
random field is a natural generalization of the familiar concept of a Markov chain, which is a 
collection of random variables with the property that, given the present, the future is (condi- 
tionally) independent of the past. If we look at the chain itself as a very simple graph and 
ignore the directionality implied by "time", then a Markov chain may alternatively be viewed 
as a chain graph of stochastic variables, where each variable is independent of all other variables 
(both future and past) given its two neighbors. A Markov random field is the same thing, only 
that rather than a chain graph, we allow the relationship between the variables to be defined 
by any graph structure, and each variable is independent of all the others given its neighbors 
in the graph. A Gibbs random field, on the other hand, is formed by a set of random vari- 
ables whose configurations obey a Gibbs distribution, which is a probability distribution that 
factorizes over all possible cliques, i.e. complete subgraphs in the graph, and the factors are 
conveniently referred to as "clique potentials". These two ways of defining a random configu- 
ration are apparently quite different pQ: A Markov random field is characterized by its local 
property (the Markovianity) whereas a Gibbs random field is specified by its global property 
(the Gibbs distribution). 

The rigorous study of the relationship between these two seemingly unrelated fields was 
initiated by Dobrushin [2] in the context of statistical physics, who considered the questions 
of existence and uniqueness of a random field subject to a Markovian conditional distribution. 
Further investigations quickly ensued. Averintsev [3] and Spitzer [4] independently proved that 
the class of two-state Markov chains is identical to the class of Gibbs ensembles on the simple 
cubic lattice. Hammersley and Clifford [5j showed that the same equivalence holds between 
a multi-state Markov field and a generalized Gibbs ensemble over an arbitrary finite graph. 
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The celebrated Hammersley-Clifibrd theorem states that each Markov field with a system of 
neighbors and the associated system of cliques is also a Gibbs field with the same system of 
cliques, and vice versa, each Gibbs field is also a Markov field with the corresponding system of 
neighbors. This implies that the joint probability and the conditional probability can specify 
each other, and serves as a theoretical basis for many modeling applications, where the global 
characteristic is captured and represented through a set of tractable local characteristics. The 
original method of proof, however, did not have great intuitive appeal, and many alternative 
proofs of this theorem were developed. Sherman [6] verified the equivalence of Markov fields and 
Gibbs ensembles under more relaxed conditions by the repeated use of the inclusion-exclusion 
principle. Preston [7] adopted a direct approach to the two-state problem and presented an 
explicit formula for the pair potential. Grimmett [8] showed that the equivalence of structure 
follows immediately from an application of the Mobius inversion theorem. A final improvement 
was done by Besag [9], who applied methods of statistical analysis and gave a much simpler, 
analytical proof of the general result. 

The nearest-neighbour Ising model in one dimension is commonly used to demonstrate the 
powerful Markov-Gibbs equivalence. Though an ordered phase only emerges at zero tempera- 
ture, this classic model is physically important in that it has a fixed point (the so-called "zero 
temperature phase transition") where the critical exponents may be sensibly defined as in higher 
dimensions. There is the astonishing empirical fact that these critical exponents depend only 
on overall features of the system, and are related to eigenvalues of the linearized renormaliza- 
tion group map near the fixed point [TO]. This universality conjecture has generated continued 
interest in the scientific community, and various approaches to the renormalization effect on 
the one-dimensional Ising model have been explored [TO] 

Consider a one-dimensional Ising model with N spins a% = ±1, labelled successively i = 
0, N — 1. We take the system size N to be very large (strictly speaking, infinite). The Gibbs 
field of this model is described by a Hamiltonian H, consisting of a nearest-neighbor interaction 
J and a magnetic field term m: 



where periodic boundary condition is imposed so that <jn = o"o ; a standard setup to ensure 
that H is translation- invariant. We focus on a specific renormalization group transformation, 
namely decimation transformation with blocking factor b. To avoid unnecessary technicalities, 
we assume that b divides N. The decimation procedure is straightforward: Fix the spins au 
for i = 0, ...,N/b — 1, and integrate out the remaining ones. This will generate a renormalized 
Gibbs field with a Hamiltonian H' having the same form as the original Hamiltonian H, but 
containing a nearest-neighbor interaction J' and a magnetic field term ml: 



The renormalized spin coefficients (J',m') and the original spin coefficients (J,m) are related 
by the decimation map: 
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6-1 6-1 
m 



= ^2 ex p I J ^ ciicji + i + ~2 ^2 ( a * + ai+i "> ) ' ( 3 ) 

iti,...,(T6_i \ i=0 i=0 / 

where C is a normalization constant. Notice that to avoid double counting, we have assigned a 
"half" of the magnetic field m (m') to each spin. 

We would like to obtain an explicit characterization of the renormalized model, but as the 
blocking factor b gets large, solving for (J 7 , m') directly from (|3J) becomes very difficult. We thus 
take an alternative approach and investigate the decimation effect on the associated Markov 
chains. As there is no finite phase transition in one dimension, we follow the common practice 
and measure the nearest-neighbor interaction strength J (J 7 ) by the Boltzmann factor k = e~ 2J 
(kf = e~ 2J ) instead. An explicit solution for (k',m') then follows from the Markov-Gibbs 
equivalence (Hammersley-Clifford theorem). The diagram below illustrates these ideas: 

original Hamiltonian ^4 renormalized Hamiltonian 

4(H) t(rv) ( 4 ) 

original Markov chain ^> renormalized Markov chain 

where: 

• (II) and (IV) indicate the Markov-Gibbs equivalence (Hammersley-Clifford theorem). 

• (I) is the decimation map on the Ising Hamiltonian (cf. ([3])). 

• (III) is the decimation map on the associated Markov chains (to be examined). 

A key tenet of the renormalization group is its explanation of universality [13]. Thus we 
would also like to verify the widely-believed universality conjecture in this special case, which 
states that the linearization of the decimation transformation with blocking factor b about 
the two-dimensional fixed point (k = m = 0) has two real eigenvalues b VT and b VH , where 
Vt = Vh = 1- Suppose we start with a Hamiltonian that is close to critical. The decimation 
map will first drive it towards the fixed point for a large number of iterations, but eventually 
will drive it away. The singular behavior of the model arises from iterating the map infinitely 
many times, and the critical properties are determined by how much time the Hamiltonian 
spends near the fixed point, when its behavior is governed by the linearization. In fact, it is 
observed that there are exact non-trivial relations between the six critical exponents (specific 
heat q, spontaneous magnetization /3, magnetic susceptibility 7, response to magnetic field at 
zero temperature 5, correlation length u, and correlation function at zero temperature rj) and 
the two eigenvalues (more precisely yx and yjj) of the linearization: 

a = 2 = 1, p = = 0, 7 = = 1, (5) 

yr yr yr 

!l " 00, u = — = l, f] = d + 2 - 2y H = 1. (6) 



d-yH yr 
(More discussions may be found in [TTJ and |14j.) 

Theorem 1.1 (Universality Conjecture). At the fixed point (k = m = 0), the Jacobian matrix 
of the renormalized spin coefficients (k',m') with respect to the original spin coefficients (k,m) 
is given by 
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The rest of this paper is organized as follows. In Section[2]we verify the universality conjecture 
by analyzing the decimation map on the Ising Hamiltonian directly (First Proof of Theorem II. ip . 
In Section [3] the statistical physics model is transformed into a probability model through the 
Markov-Gibbs equivalence (Theorems 13. II and 13. 2p . We investigate the decimation effect on the 
associated Markov chains and give an explicit characterization of the renormalized Hamiltonian 
(Theorem [33]). An alternative proof of the universality conjecture from this point of view is also 
provided (Second Proof of Theorem ll.il) . Finally, Section |4] is devoted to concluding remarks. 

2. Renormalization group approach 

In this section we will examine the renormalization group equation ([3]) directly. Although it 
is difficult to find an explicit solution to ([3]) for a large blocking factor b, the Jacobian matrix 
of partial derivatives ([7]) may be computed via implicit differentiation. 

First Proof of Theorem ll.il The decimation map ([3]) consists of 4 equations. 
1. Corresponding to o"o = °~b = 1: 

exp (C + J' + m') 



exp I J I oi + y~] (JiO-i+i + a b -i +m)](Ji+m 

,...,cr b _ 1 \ \ i=l / i=l / 



(8) 



2. Corresponding to o"o = °b = — 1: 

exp (C + J' — m') 

/ / 6-2 \ 6-1 \ 

= ^2 exp J -o\ + ^2 ViVi+i - °b-\ 1 + m y~] Oj - m 1 . (9) 

&l,—,<7b-l V V i=l / i=l / 

3. Corresponding to o"o = 1 ; °"6 = — 1: 

exp (C - J') 

= ^2 exp I J I (7i + y~] cTjO-j+i - o- 6 _i + m ^ o-j . (10) 

(T1,...,(T 6 _1 \ \ j = l / 1 = 1 / 

4. Corresponding to o"o = — 1 5 °"6 = 1 ; 

exp (C - J') 

/ / b-2 \ 6-1 \ 

= ^ exp J -ff! + ^ (jj^j + (jn + m ^ ffj . (11) 

<r 1 ,...,a b _ 1 \ \ i=l / i=l / 

Due to symmetry, (|10p and (|lip are equivalent. We may therefore assume that ([3]) breaks 
down into 3 equations: ([8]), ©, and (jlOp . To compute the Jacobian matrix of the decimation 
transformation at the fixed point (k = m = 0), we perform implicit differentiation on these 
equations at (J = oo, m = 0). As an example, we differentiate both sides of ([8]) with respect 
to m, which gives 

dC dJ' dm' . 
dm dm dm 
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exp (,J (a x + ]T (Ti(Ti+i + Im^dij f ]T a-\ 

ai,...,o b _ 1 \ \ i=l / i=l / \i=l / 



/ / 6-2 \ 6-1 

^2 exp I J ( (Ji + crjCTj+i + C7fe_i +m^(r, 

CTl,...,CT f) _ 1 \ \ 1 = 1 / 1 = 1 



(12) 



Because lower order terms become insignificant at "J = oo", it suffices to keep track of the 
"dominating terms". We have 

dC dJ' dm' 
dm dm dm 

s ( Ei=l 0"i| max ( a l + Ya=1 a i a i+l + a b-l) ) 

= 7 b2 S ' (13) 

n [vi + J2i=i o-jCTj+l + o-6-l ) 

where n(f) counts the number of a configurations that maximize f(cr), and s(5|max(/)) (which 
we will abbreviate by s(g)) is the sum of g(a) over the maximizers of f(c). Repeating this 
"dominating" procedure provides us with 6 independent equations for the partial derivatives: 

dC dJ' dm' _ 5(0-1 + ^=1^^+1 + ^-1) 

dJ + dJ + dJ ~ n{ai + ^tl w + i + ^-1) ' 

dC , Of | dm' _ s(E U «*) 1 x (15) 

dm dm dm n(a x + £ti <Wh + ^6-1) 

dC dJ' dm' = s(-oi + ^i=i QW+i - 0&-1) . . 

^ 9J ~ <-a x + Et? *W+i - a 6 _0 ' 1 j 



C7,; 



, ft*' _ , (l7) 

dm dm dm n (- ai + a^+i - <r 6 _i) 



<9C _ <9J^ _ 5((7i + X^ =1 OjOj+l - Q-6-l) 

<9J dJ n ( ai + X)»=i _ Ofc-i) 



(18) 
(19) 



It is quite clear that max(oi + <J i a i+l + °V-l) = & is achieved only when <j\ = ■ ■ ■ = 

<Jb-i = 1, and that max(— c\ +Yli=i ~ °V-i) = bis achieved only when a\ = ■ ■ ■ = <Jb-i 

1 . The harder task it to determine when max(<Ti + E^~^ 0^0^+1 — C6-1) is obtained. Because of 
the single "— " sign, it is not possible for all the b terms in this sum (01, 01 02, 0-6-2 c&-i, — Cb-i) 
to be 1 as in the previous two cases. An ideal maximizer should have b — 1 terms with value 1 
and only one term with value —1. We claim that each one of the b locations of —1 corresponds 
to exactly one a configuration: Suppose the ith term has value —1 (01 = — 1 for i = 1, 
Oj-iOi = —1 for 2 < % < b — 1, or cr^-i = 1 for i = b), then we must have a\ = ■ ■ ■ = o^_i = 1 
and Oj = • • • = db-i = —1- (fl4j) — (fT9j) are thus simplified: 



d(7 + ar + cW_ & ^L + ^ + ^L- b (20) 

5J <9J <9J 9m 5m 5m 

5C ^_^__ b ^_ 9J^_dm^ _ b (21) 

dJ dJ dJ dm dm dm 

— - — -b-2 —- — -0 (22) 

dJ dJ dm dm 
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Solving USED— H22D yields 

— = 1 — = — = — = b (23) 

dJ dm dJ dm 

which further implies that the Jacobian matrix Jac ([7]) is diagonal, i.e., = A^- = 0. To 

complete the proof of the universality conjecture, it remains to verify that = b. We perform 

the "dominating" procedure as before. For notational convenience, we temporarily denote 

°"i + Si=i <W+i - CTft-l by f(a), and o x + YaZI crm+i + 0"5-i by 5(0-). Dividing ([in]) by © 
at the fixed point ( J = 00, m = 0), we have 

k> = exp(-2J') = ^(/)exp(J-max(/)) 
n(<?) exp (J • max(<7j) 

= t «p((t-2)J) =t eM _ 2J) = bL (24) 
exp(oJ) 

□ 

3. Markov chain approach 

In this section we will transform the statistical physics model into a probability model and 
investigate the decimation effect on the associated Markov chains. This is a special case of 
Hammersley-Clifford theorem where the exact correspondence between the Markov field and 
the Gibbs field may be worked out explicitly. The idea is to regard the Ising system as a 
two-state Markov chain with transition probability matrix 

"-(^'i-J- (25) 

where 

p = P(ffi = l|<7 = -l), (26) 

and 

q = P(<7i = -l|a = 1). (27) 

Theorem 3.1 (Hammersley-Clifford). The Ising Hamiltonian H |IJ) is fully characterized by 
the transition probabilities p and q. 

Remark. The transition probabilities p and q and the spin coefficients k and m are related by 
t36\), (37), |^i[ ), and The spin coefficients fixed point (k = m = 0) thus corresponds to 



the transition probabilities fixed point (p = q = 0). 

Proof. Baxter [TJ] showed that the mean and covariance of the Ising spins in the infinite- volume 
limit are functions of the spin coefficients k and m: 

„ sinhm . , 

Eao = (28) 

V sinh m + k 2 



k 2 cosh m — v sinh 2 m + k 2 

Cov(o- ,o-i) = — - 2 — =. (29) 

smn m + K cosh m + y s i n b m + k 2 



Through the Markov-Gibbs equivalence, we show that (I28h and (|29() may alternatively be 
viewed as functions of the transition probabilities p and q. Recall that the Markov chain has a 
stationary distribution: 

P(a = -1) = P(o- = 1) = (30) 

p+q p+q 
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which readily displays the dependence of the mean on the transition probabilities, 

Ea = Eat = — • (31) 
P + q 

To obtain an analogous expression for the covariance, we resort to the tower property of con- 
ditional expectation, 

Ecr CTi = E(c7 E(<Ti|cro)) 
= P(<ti = l\a = l)P(a = 1) " P(ffi = -Iko = l)P(ffo = 1) 
-P((7i = l|ff = -l)P(ffo = -1) + Pfa = -lko = -l)P(o- = -1) 

(P-^) 2 . / n \ 4 P<? , qo x 

= +{1 -"-' ,) JFTW (32) 

which then gives 

4pq 

(p + q) 

The two characterizations of the Ising Hamiltonian H (JTJ) are thus connected by: 

p — q sinh m 



Cov(cjo, oi) = Ea ai - Ecr E<7i = (1 — p — q) y- ■ ^ 2 . (33) 



P + q Vsinh 2 m + k 2 ' 



(34) 



^ 4pg A; 2 cosh m — y sinh 2 m + fc 2 . 

(P + <?) 2 sinh 2 m + & 2 cos h m + Vsinh 2 m + k 2 



It is not hard to derive an explicit expression of p and q in terms of k and m from (|34p and 
351): 



V sinh 2 m + k 2 -\- sinh m „ , 

P = ~ , 2 (36) 

cosh m + v sinh m + fc 2 



\/ sinh 2 m + k 2 — sinh m , _. 

q = 1 2 (37) 

cosh m + V sm h m + k 1 

The reverse direction, however, requires more work. For computational convenience, we make 

a change of variables, A = sinhm, B = sinh 2 m + k 2 . Then f)36[) and (j37j) become 



Vb + A = pv / A 2 + 1+ P Vb, (38) 



V^B - A = q\/A 2 + \ + qVB. (39) 

Dividing (|59"j) into (|38|). we have 

-yl - gv^ 2 + 1 9-1" 

This is an equation for ^4 only, and an explicit expression of k and m in terms of p and q follows 
easily: 

m = I log ( : J— ^ . (42) 

□ 
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Theorem 3.2. The renormalized Ising Hamiltonian H' (dj) is fully characterized by the renor- 
malized transition probabilities p' and q' , where 

p> = P(a b = l\a = -1), (43) 

and 

q' = P(a b = -l\a = 1). (44) 

Remark. The renormalized transition probabilities p' and q' and the renormalized spin coeffi- 
cients k' and m! are similarly related as in A36\), ft37\ ), ft41\ ), and ft42\ ). 

Proof. This follows from Theorem 13 . 1 1 once we realize that site and site b are nearest neighbors 
after decimation transformation with blocking factor b. The 6-step transition probability matrix 
P b represents the decimation map on the associated Markov chains, and is given by 

b 



(45) 








( p i 


vl ? ] ] 








p+q J J 



1 -1 

p+q p+q 



where the first equality is simply the spectral decomposition of the matrix P. This then implies 
that 

p'=^-(l-(l-p-q) b ), (46) 
p + q 

and 

q' = ^-{l-{l-p-q) b ). (47) 
p + q 

□ 

Theorem 3.3. The decimation map ([3J) identifies the connection between the renormalized 
Hamiltonian H' ^) and the original Hamiltonian H ([!]). 

Proof. We follow (II), (III), and (IV) as shown in @. The original Ising model is described 
by a Hamiltonian H with spin coefficients k and m. (II) indicates the alternative view of this 
system as a two-state Markov chain with transition probabilities p and q (cf. (|36p and (|37[) ). 
(Ill) then transforms this Markov chain into a renormalized Markov chain with renormalized 
transition probabilities p' and q' (cf. (|46p and (|47p ). Finally, (IV) recovers the renormalized 
spin coefficients k' and m! of the renormalized Hamiltonian H' (cf. (|4ip and ([42])). 

(fc.m) ^ (p,q) ( ™ } (p',q') ^ (k',m') (48) 

□ 

Second Proof of Theorem \l.l[ Theorem l3.3l establishes an explicit expression of the renormalized 
spin coefficients k' and m! in terms of the original spin coefficients k and m (cf. (|48p ). To 
evaluate the Jacobian matrix Jac (JT]) at the fixed point {k = m = 0), we start by considering 
^jr- and with m held fixed at zero. By (II), on the m = curve, 

" = " = TTk- < 49 » 



(III) then gives 

n' = a' = 

2 



p' = q' = \[l-(l-2p) b ), (50) 
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which further implies, by (IV), that 

« = i^. ««> 

rri = 0. (52) 
We conclude that = from (|52j) . and by applying the chain rule to (Tl9j) . (j50j) . and (|5Tj) . 
that ^ = 6. 

We proceed with the calculations for and |^ with /c held fixed at zero. By (II), on the 
k = curve, either p or q is zero, depending on the sign of m. Without loss of generality, 
assume m > 0. In this case, 

2 sinh m , . 

P= 7—r -r-r— , (53) 
cosh m + sinn m 

g = 0. (54) 



(III) then gives 

which further implies, by (IV), that 



p' = l-(l-p) b , (55) 
q' = 0, (56) 

k' = 0, (57) 

m' = -Uog(l-p'). (58) 

We conclude that |^ = from ([57|) . and by applying the chain rule to ([53]) . ([55]) . and ([58]) . 
that £ = 6. □ 

4. Concluding remarks 

This paper aims at an explicit characterization of the renormalized Hamiltonian after dec- 
imation transformation of a one-dimensional Ising-type Hamiltonian with a nearest-neighbor 
interaction and a magnetic field term. We transform the statistical physics model into a prob- 
ability model through the Markov-Gibbs equivalence and analyze the decimation effect on the 
associated Markov chains. As the Ising model is a prototype for a wide variety of spin models, 
it is expected that the exploitation of Markov-Gibbs equivalence in this special case will shed 
light on the application of renormalization group ideas in a more general setting. Two different 
proofs of the universality conjecture are presented, one based directly upon the renormalization 
group equation, and the other from the Markov chain point of view. Although the first proof 
does not employ advanced mathematical methods, it provides a new perspective on the renor- 
malization flow. For example, it has been verified, following similar ideas, that one-dimensional 
(/-state Potts model (q > 2) exhibits the same eigenvalue statistics yx and yn, independent 
of the number of states q and the blocking factor b. (The percolation limit q — > 1, however, 
remains open, and is believed to display different critical features.) The second proof uses ideas 
from Markov chains, and is expected to work with higher-dimensional g-state Potts models as 
well, where the covariant matrices may be expressed in terms of the random cluster represen- 
tation of Fortuin and Kasteleyn |15| . As the number of dimensions d and the number of states 
q get large, it will be harder to write down exact formulas for the transition probabilities in 
the covariant matrices, but the Metropolis and Glauber algorithms should provide a reasonable 
approximation scheme. Since Markov chains may take both discrete and continuous values, an 
advantage of exploring this second perspective is that we can also consider decimation with 
spin scaling applied to continuous spin systems, not just discrete systems like Potts and Ising 
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models, and hence avoid the lack of spin rescaling, a common problem encountered in a "pure" 
decimation. In summary, we hope this rigorous investigation will provide insight into the in- 
trinsic structure of the renormalization group transformation and help us better understand 
the nature of universality. 
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